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Abstract 


The  Helmholtz  integral  equation  may  be  formulated  by  combining  the  scalar  wave 
afiution  and  Euler’s  equation  for  motion  unthin  a  fluid.  The  solution  to  this  integral 
equation  yields  the  radiated  pressure  from  a  submerged  vibrating  body.  The  solution 
obtained  pom  applying  Hamilton's  variational  principle  governs  the  response  to  a  given 
surface  pressure  e'xitation.  These  two  solutions  may  be  used  in  conjunction  to 
investigate  the  scattering  of  incident  sound  waves  pom  shell-like  bodies.  In  this  report, 
the  scattering  pom  a  submerged  steel  spherical  shell  having  a  thickness  to  radius  ratio 
of  0.01  is  investigated  for  ka  values  up  to  S. 
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Sound  Wave  Scattering  from  an  Elastic 
Spherical  Shell 


1.  Introduction 

Anechoic  coatings  have  the  abUity  to  modify  the  scattering  behaviour  of  sound 
waves  incident  upon  a  body.  The  first  step  in  the  investigation  of  anechoic 
coatings  is  to  consider  their  effects  when  applied  to  simple  objects  such  as  a 
sphere  or  a  spherical  sheU.  Beft»«  an  analysis  of  this  is  possible,  it  is  necessary 
to  have  a  firm  grasp  of  the  scattering  characteristics  of  the  body  prior  to  the 
application  of  the  coating. 

The  purpose  of  this  report  is  to  extend  the  analysis  given  in  Partridge  {11  for 
the  rigid  sphere,  to  the  scattering  behaviour  of  a  thin  elastic  spherical  shell. 

The  mathematical  derivation  of  the  scattered  pressure  requires  the  use  of  both 
the  Helmholtz  integral  equaticm  and  Hamilton's  variational  principle.  This 
attalysis  essentially  follows  that  given  in  the  work  of  lunger  and  Feit  [2]. 
Combining  this  with  various  scattering  parameters  in  the  literature  [3*71,  such 
as  target  strengfo,  back>8cattering  cross-section,  reflection  factor  and  intensity 
field  around  the  body,  Uie  overall  behaviour  of  the  elastic  shell,  nuy  be 
analysed.  In  a  subsequent  report,  these  results  will  be  compared  with  those 
obtairted  from  coated  spherical  shells. 


1.1  Helmholtz  Integral  Equation 

A  sourvl  wave  can  be  thought  of  as  a  tinne  dependent  pressure  fluctuation  P 
around  the  static  pressure  P,  in  a  compressible  fluid  such  as  water.  The 
temporal  and  spatial  variations  of  this  fluctuation  ate  governed  by  the  wave 
equation 


where  c  is  the  wave  propagation  speed  through  the  fluid.  If  we  restrict 
ourselves  to  steady-state  situations,  where  the  pressure  varies  periodically  with 
time,  i.e. 


P(i.t)  -  P(r)e' 


then  the  above  reduces  to  the  "Helmholtz"  or  "steady-state"  wave  equation 


(V^  ♦  k^)P{r)  -  0 


where  dte  acoustic  wavenumber  is  given  by  <:  =  uj/c.  Euler's  equation  of 
motion  within  a  non-viscous,  homogeneous  fluid  is  given  by 

' -P 


where  p  and  v  are  the  density  and  velocity  respectively.  For  a  non-viscous 
fluid,  pressure  can  act  only  in  a  direction  normal  to  the  surface  boundary 
and  letting  the  normal  component  of  velocity  to  the  boundary  be  given  by 


where  to  deflnes  the  surface  the  above  boundary  condition  reduces  to 


■  (iuplil;  on  S 


where  n  is  the  unit  outward  normal  to  S^.  The  pressure  radiated  from  a  body 
vibrating  with  a  surface  velocity  distribution  w,  is  governed  by  equations  (1 J) 
and  (1.6).  From  Partridge  (IJ,  these  two  equations  rr\ay  be  combined  to  form 
the  Helndtoltz  Integral  equation 


r  outsides^ 


where  Green's  function  G(r  |  f„)  represents  the  field  at  a  point  r  due  to  a  unit 
impulse  located  at  r^.  If  the  Green's  function  satisfies  the  Neumarm  boui  •  iaiy 
condition 


00 

^  =  0  on  the  boundary 


then  the  Helmholtz  integral  equation  reduces  to  the  simple  integral 
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Pit)  -  (iwp)  //[G(£jr„>t^(,„)]  <fS„ 


1.2  Pressure  Field  for  a  Spherical  Radiator 


For  a  vibrating  spherical  body  of  radius  a,  the  above  integral  becomes 


F(r,d)  •  (uopUlKa^)  sind„  d-dg 


The  Green's  function  and  surface  velocity  are  given  (1]  by 


””  hi,  ity) 

G(r,«|fl,«^)  •  — L_  £  (2n  ♦l)P„(cosd)P„(cosd„)  _  (1.11) 

4itfaj^  n«o  h'(ka) 


K(ka) 


ii;(d)  -  £  W„P„(cosi3) 


where  P^  (cos  9)  and  /i„  (br)  are  the  Legendre  and  Hankel  functions  respectively 
and 


^  (2n^l) 

*V_  ■ _ 

n  2 


J  P^fnlwCTlldil,  Ti  »  cos^ 


Substituting  equations  (l.ll)  and  (1.12)  into  (1.10),  yields  the  pressure  radiated 
from  a  spherical  body 

“  ,  hAkr) 

P(r,d)  -  i(pc)  £  W„P„(cos6) -I -  (l.h 

n*o  h„  (ka) 


This  solution  will  be  applied  later  to  obtain  the  scattered  pressure  from  a 
spherical  shell  impact^  plane  waves. 


1.3  Radiation  Loading 

When  a  submerged  body  vibrates,  it  sets  up  pressure  fluctuations  in  the 
surrounding  fluid.  These  fluctuations  (termed  the  radiation  loading)  then  act 
back  on  the  structure.  This  radiation  loading  modifies  the  forces  acting  on  the 


I 
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structure,  and  since  these  pressures  depend  on  the  velocity  at  which  the 
structure  vibrates,  a  feedback  coupling  between  the  fluid  and  the  structure 
exists.  The  ratio  of  surfoce  pressure  to  the  surface  velocity  is  defined  as  the 
specific  acoustic  impedance.  In  modal  form  this  becomes 

z  •  f  /iV  <5 

‘•n  In'  "n 

where 


5^  /„P„(COSd),  taffi)  «  53  W„Pn(cOS-d) 
n«o  n-o 


(1.16) 


The  specific  acoustic  impedance  measured  on  the  surface  of  the  source, 
determines  the  radiation  loading,  and  may  be  expanded  into  its  real  and 
imaginary  components  which  are  designated  by  and  (om^.  Therefore, 
equation  (1.15)  may  be  rewritten  as 


«  r_  W_  ♦  m„  IV. 
n  n  n  n 


(1.17) 


The  surface  pressure  is  the  sum  of  two  terms,  (Da  term  proportional  to  the 
surface  particle  velocity  w  embodying  a  resistive  or  damping  force,  which 
represents  the  energy  lost  by  the  pulsating  sphere  in  the  form  of  radiated 
acoustic  energy,  and  (2)  a  term  proportional  to  the  surface  particle  acceleration 
associated  with  the  mass  of  fluid  set  into  motion  by  the  pulsating  surface  of  the 
spherical  source.  Using  equation  (1.14)  for  the  radiated  pressure,  the  modal 
acoustic  impedance  for  a  spherical  body  becomes 


z 


n 


f(pc) 


h^ika) 


rn-itom„ 


(1.18) 


where 


r 


n 


(pc)Re 

(pc)l(*B)|)i^(to)|l-2 

(pcXto)^"^ 

(rfl)^ll.3...(2n-l)l2' 

(pc). 


(Jbj)^  «  |2n-l  I 
(far)  >>  «^  +  l 


(1.19) 


and 


10 


2.  Equations  of  Motion  for  a  Spherical  Shell 

2,1  Formulation  of  Problem 

The  solution  to  the  equations  of  motion  for  an  air-iilled  spherical  shell  subjected 
to  an  excitation  pressure,  yields  the  velocity  distribution  around  the  shell's 
surface.  The  surface  pressure  exdtes  the  normal  vibration  modes  of  the  shell. 
The  exterior  boundary  pressure  and  the  displacements  on  a  spherical  surface 
(called  the  shell  middle  surface)  which  is  located  centrally  between  the  physical 
surface  boundaries  of  the  shell,  may  be  expressed  as 

«• 

P,(.a,-d)  *  2}  /„P„(T1),  n  -  cosd 

n«o  (2.1) 

m  m 

uW  *  23  li„(l dP„(Ti)/«/Ti,  u;(d)  -  £ 

n"o  n»o 

where  u  and  zu  are  the  "tangentUl”  and  "nomnar  displacement  components 
respectively.  The  velocity  distribution  is  related  to  the  displacements  via 
u  =  (-  I'colu  and  w  =  (-  iQ>)ur.  In  order  to  derive  the  equations  of  motion,  several 
assumptions  must  be  made: 

(a)  Displacements  are  assumed  small  compared  vrith  the  shell  thickness. 

(b)  The  shell  thickitess  is  small  compared  with  the  radius  of  curvature  of  the 
sphere. 

(c)  Fibres  of  the  shell  initially  perpendicular  to  the  shell  middle  surface 
remain  so  after  deformation  and  are  themselves  not  subject  to  elongation. 

(d)  The  normal  stress  a„  acting  on  planes  parallel  to  the  shell  middle  surface 
is  negligible  in  comparison  %vith  other  stresses. 

These  assumptions  allow  us  to  reduce  the  problem  of  shell  deformation,  to  the 
study  of  middle  surface  deflection.  The  equations  of  motion  are  derived  from 
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displacements.  Hamilton's  principle  states  that  if  the  work  done  on  the  shell 
by  the  extenud  pressure  causes  the  displacemonts  to  undergo  a  variation 
between  the  times  and  fj,  then 


(2.2) 


where  £j.,  £,  are  the  kinetic  and  strain  energies,  and  5a;  3  the  work  done  by  the 
pressiue.  Consider  a  thin  spherical  shell  of  thickness  h  and  mean  radius  "a"  as 
shown  in  Figure  1.  In  this  derivation,  only  non-torsional  axisymmetric  motiorts 
of  the  shell  ate  considered,  making  the  variables  independent  of  the  angular 
coordinate  0.  The  shell  displacements  (i,  and  W,  can  be  written  in  terms  of 
the  middle-surface  displacements  u,  and  w  using  {2] 


U. 


1  *1 


f  ^ 

X 


9u) 


W,  =  a; 


(23) 


where  x  =  r  -  a.  The  total  kinetic  energy  of  the  system  may  be  written  as 


*  jPs Iff 

2it  K  ♦h/2 

=  ip,  /  /  f  (iJ,^  +  IV3  )  Sind  ifd  rfi  (2.4) 

^  00  -in 

» 

=  KQ^ha^  j (li^  +  w^)  Sind  id 


In  spherical  coordirtates,  the  strains  and  stresses  are 


1 


(a*x) 


dU, 

IF 


■W, 


du 


+  w 


f  da  d^a; ) 


(25) 


i _ _  (cotdl/j  +  W.)  »  1.  (cotdu  +w)  *JL  cotd 

a 


(a*x)  s  -s  a  '  .2 


a  - 


9w 

Id 


- 5- 

d-’/) 


(1  -v^) 
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Scafiered 
spherical  wave 


Incident 
plane  wave 


Figure  1:  Sound  sattering  from  an  elastic  spherical  shell 


The  strain  energy  of  the  system  in  terms  of  these  components  is 


dV 


^  *h/i 

2  /  /  J 


o  0  -h/2 

1C 


Jt£h 


/ 


(1-V^)  o 
_ ^  " 


'•  'll 

^  +11?  ♦  (coti&u  *w^  * 


2v 


^u 


^  tcEhfi^  1’ 

/ 1  _  .»2\  J 


*w 

\  ~  J 


37 


+  2vcoti) 


f  3u  _  3^x0 


(  \ 

du? 


u  - 


da? 


■55 


icottu  *w) 


sill’d  d-5 


V 


■55 


sin-6  did 


(2.7) 
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where  the  contribution  of  the  bending  stresses  to  the  strain  energy  is 
represented  by  the  term  ■  k^/\2a^.  The  work  done  by  tfte  surface  pressure 
(applied  normal  to  the  uttdefonned  shell  surface)  is 


-  J/p.wdS  -  J  Jp.  w  Sind  da 

O  O 

K 

■  (2jifl^)  JPjjU;  Sind  dd 


Applying  Hamilton's  principle  with  respect  to  the  variations  du  and  6v  to  obtain 
the  eqtrations  of  motion  for  the  spherical  shell  yields 


d^u  du 


(1  ♦p  )  — -  >cotd  -(v+cord)M  -  P 


I  r,2 


-P^cotd  »(<!  *v)  »p^(v  ♦cot^d))  -  f  ^  * 


■5d  ~r 


p^  i_4  +2p^cotd  *  P^cot^i») 

dd^  d^ 

♦cotd  ((2-vcot^d)p2  -(1  ♦vllw-p^  i!!f!-2p^cotd  ^ 

dd^  dd^ 

+  P^(l  +  v  +  cot^d)  -p^cotdil-v  +  cot^d)^  -2(1  ♦v)u)  -  “  f  - 


-P,(l-v^)fl^ 


w*'ere  the  wave  speed  in  the  shell  material  Cp  is  given  by 


p^d-v^) 

Introducing  the  variable  rj  =  cos  9  and  assuming  harmonic  time  variation,  the 
foregoing  set  of  equations  may  be  rewritten  in  the  form 


f.uu(«)  ♦  *  0 

L^(«)  .  L^(a.)  .  n^u-  -  -P. 


where 
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(1  (1  -Tl^ (1  -V) 

[  ^ 

^  f  ^  \ 

(l-TJ^)’''^  I3^(l-V)-(1  +V)l4- 

I  In  7?  ’’J 

-  ip^d  -v)-(i  ♦v)i^{i  -n^)’''^*p^v?4-d 
^  dn  on 

-p^vl  -p^(l-v)V^  -2(1  ♦v) 


Substituting  u  and  w  for  their  respective  series  representation,  yields 


where  A,  =  n(n  +  1).  Therefore  to  obtain  the  response  of  the  shell,  the  above 
linear  equations  are  solved  for  each  n,  to  obtain  U„  and  W„  and  then 
substituted  into  equations  (2.1)  for  u  and  to  to  yield  tlw  middle  surface 
displacements  of  the  shdl  due  to  the  exdtation  pressure  P,. 


2.2  Natural  Frequencies  of  the  Shell 
2.2.1  In  a  Vacuum 

The  resonant  frequencies  for  the  free  vibration  of  a  spherical  shell  may  be 
obtained  by  setting  /„  to  zero  in  equation  (2.14),  and  then  eliminating  L/^  and 
W^.  This  yields  the  frequency  equation 

a*  -  (d  OvA^)  -p^d  -v-A^  -vA^))tf  dA,  -2)d  -v')  ^2.15) 

♦  p2(A^-4A^+Aj,(5-v2)-2(l-v2))  -  0 

where  the  non-dimensional  frequency  Q  is  defined  as  (lea)  (c/Cp)  and  Cp  is  given 
by  equation  (2.10).  The  above  equation  may  be  written  in  the  form 
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-  (£4”)^j  [n^  -  (o^V]  -  0 


(2.16) 


where  0„®  are  the  two  resonant  frequence,  of  which  Q„®  is  the  huger. 
For  e«:h  mode  chancteiued  by  a  specific  value  of  n  >  0,  there  are  two  distinct 
positive  real  resonant  frequencies,  but  for  n  s  0  correspcmding  to  what  is  called 
the  breathing  mode,  there  is  <miy  one  real  positive  foe^ency,  the  other  beii^ 
imaginary. 


2.22  Submerged  in  Viater 

In  order  to  determine  the  resonant  frequencies  of  the  submerged  shell,  we  treed 
to  go  back  to  the  equations  (2.14),  i.e. 


lrf-(l +p2)(v+X^-l)JLr„  ♦V)]W„  -  0 


(2.17) 


-\,{p2(v»A^-l)>(l>v)JLf„*[rf-2(l»v)-p2(vX„-l))W„  -  - 


^(l-v^) 


TT 


To  determine  these  frequencies,  a  radiation  damping  term  must  be  introduced. 
The  loading,  in  the  absettce  of  a  driving  force,  is  expressed  in  terms  of  the 
spediic  acoustic  impedance  as 


P.  -  -  E  "„^V„P„(cosf)) 

^n-o  (2.18) 

*  5^  [(iB>)r„  ♦(u>^)m„J  W„  Pn(co5«) 

n«o 


The  negative  sign  arises  because  the  applied  pressure  was  defined  previously  as 
positive  outward,  whibt  the  radiation  loading  is  positive  inward.  Therefore, 
replacing  /„  in  equation  (2.17)  with 

-  -l(ru))r„  ♦  ((B^)m^)  W„  W.19) 

yields 


(Q^-(l+P^)(v*\,-l)HJ„-ip2(v*\,-l)*(l*v)lW„  -  0  (2.20) 


!♦ 


kini 


-2(l*v)-p^(v+k„-^> 


|W„«0 
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Combining  equaaons  (2^0),  and  noting  that  the  natural  frequencies  are  dte 
roots  of  the  real  compoitent  of  the  resultant  equation  yields 


a^O*  *ajCi^  *02  •  0 


(2^1) 


where 


tfo  -  n  > 
a,  - 

Oj  -  [2(l+v) 

The  solution  to  this  equation  yields  the  natural  frequencies  for  a  submerged 
shell.  It  is  a  transcendental  equation  since  m„  given  by  equation  (1.20)  is 
frequency  dependent. 

Listed  below  are  the  values  that  %vill  be  used  in  this  report  for  both  the 
propagation  medium  and  tl^  spherical  shell. 


Propagation  medium 


Sound  velocity  in  water 

c 

= 

1500 

(m/s) 

Water  density 

P 

S 

1000 

(kg/m^) 

Spherical  shell  (steel) 

Mean  radius 

a 

1,0 

(m) 

Shell  thickness 

h 

= 

0.01 

(m) 

Density 

Ps 

= 

7700 

(kg/m^) 

Young's  modulus 

E 

= 

195  X  10'° 

(Pa) 

Poisson's  ratio 

v 

X 

029 

(•) 

The  rtatural  frequencies  for  the  above  steel  shell  in  a  vacuum  and  submerged  in 
water,  are  tabulated  in  Ibble  1.  From  the  table,  it  may  be  seen  that  radiation 
loading  reduces  the  rratural  frequencies  of  the  shell. 
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Table  1:  Effect  of  submergence  on  the  normal  modes  of  a  spherical  shell 


Branch 

n 

Dimenaionless 

Natural  Frequency  Ci 

In  vacuo  Submerged 

Dimensionless 

Frequency  parameter  ka 

in  vacuo  Submerged 

Lower 

nn‘" 

0 

imaginary 

imaginary 

0.00 

0.00 

1 

ODOO 

0.000 

0.00 

0.00 

2 

0.705 

0225 

2.47 

1.14 

3 

OSU 

0.411 

2.92 

1.44 

4 

0S86 

0.468 

3.10 

1.64 

5 

0.912 

0314 

3.19 

1.80 

6 

0.930 

0348 

3.25 

1.92 

7 

0.944 

0582 

3.31 

2.04 

8 

0.958 

0.610 

326 

2.14 

9 

0.975 

0.639 

3.41 

2.24 

Upper 

fln® 

0 

1.606 

1.156 

1 

1.967 

1.763 

2 

2.715 

2267 

3 

3260 

3.417 

4 

4392 

4.444 

5 

5372 

5.460 

6 

6359 

6.471 

7 

7350 

7.478 

8 

8344 

8.483 

9 

9339 

9.487 

2.3  Velocity  Response  of  Shell  due  to  an  Arbitrary 
Surface  Excitation 

23.1  In  a  Vacuum 

The  middle-surface  displacements  caused  by  a  surface  excitation 

mm 

P,W  *  /nP„(COSlS)  (222) 

n«o 

were  earlier  expressed  as 

mm  mm 

«(«)  -  5}  U„(l (fP„(Tl)/(lTl,  wfd)  -  52 

nao  n«o 
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where  the  coefficients  of  li„  and  are  given  by 

(2^4) 


Solving  the  above  for  I4  and  yields 


„  [a*(l-v^L  [PV*\-i)*a*v)l 

•’I  flk  r  (0»-D‘«)*(Q»-Q.»)> 


a*(l-v*)L  rQ^-(l*P^(v*;i.-l)l 

^  1  “  (Q*-Qi'V(0*-0f)* 


The  above  may  be  simplified  somewhat  by  introducing  a  term  known  as  the 
modal  medumical  imp^ance  Z„,  which  is  defined  as  the  ratio  of  the  resultant 
surface  pressure  to  the  normal  velocity,  i^. 

Z„./„/{-ia)W„) 

.fp,Cph]in?-(Qj,*Viirf-(n^Vl  (22< 

J  tQ^-(l+p^)(v*X  -1)1 


Then  for  a  given  excitation  pressure  the  normal  velocity  may  be  expressed  as 


•  r  (/n/Zn>Pn<C'»^> 


2.32  Submerged  in  Water 

Let  the  excitation  and  resultant  pressures  be  represented  by 


-  r  "  E  /nPn(«>S^> 


Then  the  resultant  pressure  may  be  written  in  modal  form  as 


/n  -  /n 
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The  resultant  pressure  is  however  also  related  to  the  modal  mechanical 
impedance  via  2„  =  Combining  these  two  expressions  yields  the 

resulting  normal  velocity,  i^. 


d;(6)  -  £ 
n*o 

where  and  are  given  by  equations  (L18)  and  (2J26). 


3.  Scattering  of  a  Plane  Wave  from  a 
Spherical  Shell 


- —  P_(cosd) 


(230) 


Consider  a  distant  point  sound  source  which  generates  a  continuous  sound 
wave.  Far  away  from  this  source,  and  over  suitably  restricted  regions,  these 
waves  may  be  said  to  approximate  plane  waves.  Consider  these  plane  waves 
to  be  incident  upon  a  stationary  thin  elastic  spherical  shell,  whose  surface  5^  is 
given  by  the  position  vector  shown  in  Figure  1.  The  following  analysis 
differs  from  that  of  a  rigid  sphere,  in  that  the  scatterer  is  allowed  to  deflect 
elastically.  In  this  situation,  Ute  structure's  dynamic  response  radiates  a 
pressure  field  which  further  modifies  the  sound  field. 

If  we  define  a  wavenumber  vector  k,  whose  magnitude  is  k  and  which  lies  in 
the  direction  of  wave  propagation,  then  the  incident  pressure  wave  may  be 
written  as 


(3.1) 


In  spherical  coordinates,  a  wave  incident  from  the  =  180*  direction  is 
represented  by 


P((r,d)  -  Pgexpidrrcosd)  (3-2) 

If  the  pressure  scattered  by  an  elastic  boundary  is  denoted  by  P^  then  the 
resultant  pressure  for  a  thin  shell  becomes 

P  -  P,*P„  (33) 

The  term  may  be  broken  up  into  two  components,  a  term  P^  which  is  the 
pressure  scattered  if  the  shdi  Irehaved  like  a  rigid  body,  and  another 
component  P,,  i.e. 
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-  P,,  ♦  P,  (34) 

Therefore,  the  resultant  fmssure  becomes 

P  -  P,  ♦  P,.  ♦  P,  (33) 

From  equation  (1.6),  the  resultant  pressure  and  velocity  are  related  by 

•  (i(Op)tl/  on  Sg  (34) 

where  the  resultant  velodty  on  Sg  is 

lb  ■  (ti»j  ♦  ti>,.)  ♦  a;,  (3.7) 


Substituting  for  w  and  P  in  the  above,  yields  the  four  relations 


w  •  tbj 

(a) 

Wj  +  »  0 

(b) 

II 

II  i. 

(0 

(d) 

In  order  to  apply  the  Helmholtz  integral  equation  solution,  the  scattered 
velocity  comporwnts  are  represented  in  the  form 


(34) 


“  IT  (cos4),  (cos«)  <^.9) 

n-o  n«o 

The  incident  wave  may  be  represented  in  tenrts  of  a  series  of  concentric 
sfrfterical  waves  as 


P,(r,4)-P,  5^  l(2n*l)(i)'’;„(]kr)]P„{cos©)  «-W> 

n«o 


where  /„  (kr)  is  the  spherical  Bessel  function.  From  equations  (34b,  d),  the 
coefficients  of  may  be  obtained,  i.e. 

K”  -  (iPo/pc)(2n  ♦])(»■)" ;'(*»)  <3- 

In  order  to  apply  the  Helmholtz  integral  equation  solution,  P^  is  regarded  as 
the  pressure  radiated  by  a  spherical  source  vibrating  with  velocity 
therefore 
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(3.12) 


„  A  1  •-  h  (kr) 

-  i<pc)  £  IV*  P„(co9tf)  J1 - 

— f  *  «  f.  _  V 


)t„  (Jta) 


A  n  / 

-p*  E  (2n^l)(i)";„(lka)P„{cos«) -2 - 


h„(ka) 


The  Helmholtz  integral  equatitm  solution  can  also  be  used  to  determine  P^ 
Equation  (3.8c)  indicates  ^t  P,  has  the  nature  of  a  radiated  pressure  caused  by 
a  sphoical  body  vibrating  with  velocity  Wj.  This  vibration  is  due  to  the  elastic 
structure  respoitding  to  the  surface  pressure  generated  by  the  incident  wave 
field.  This  surface  pressure  may  be  written  as 

P,(d)  -  P,(a,d) +  P,.(a,d) 

T'  <3.13) 

-  E  /n  Pn<‘^os^> 


where  from  equation  (3.10)  and  (3.12) 


P„(2n  ♦1)(J)"' 
(kai^h^  ika) 


From  equation  (2.30)  in  the  previous  chapter,  the  modal  velocity  W„  caused  by 
the  modal  excitation  /„*  is 


-  -Po  E 


(2n  +  l)(i)"*’ 


«-o  (z„  ♦Z„)(fai)^lhn(te)l 


P„(cosO) 


From  the  Helmholtz  solution,  the  radiated  pressure  P,  is 

P,(r,d)  -  i(pc)  E  Pn(cosO)  . . . . 

n-v  i^(ka) 


Substituting  for  yidds 
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P,(r,i9)  -  (pc)P^  £  - !1 -  P„(co»d)  (3.18) 

n-o  (2„*Z„)I(lbj)h^(te)]2 

When  r  is  a  long  distance  from  the  boundary  i.e.  "far-fidd",  U«  Hankel 
function  Oar)  may  be  approximated  by 


exp(tiT)  (-t)"** 


and  the  scattered  pressure  fields  become 


expi(fcr) 
kr  D«o 


(2n*l)P„(cosiJ);'(fcfl) 

- P„  (costs) 

hl^ika) 


(3.19) 


P,(r,«) 


-iipOP^ 
- 1? - 


expi(jtr) 


E 


(2«  +  l)P„(cos«) 
(z„-Z„)l(lai)fc'(fca)l2 


(3.20) 


4.  Methods  of  Presenting  Scattering  Solutions 


To  examine  the  scattering  behaviour  of  a  thin  elastic  spherical  shdi,  certain 
parameters  for  the  scattering  body  need  to  be  determined.  These  parameters 
indude  (i)  the  resultant  pressure  field  in  the  backscatter  direction,  (ii)  the 
backscatteiing  cross-section  of  the  target,  (iii)  the  target  strength  of  the  ^tterer, 
(iv)  the  reflection  factor  around  the  scattering  body  and  (v)  the  intensity  held 
around  the  scattered  To  determine  these  parameters,  values  previously  given 
in  chapter  2  for  the  properties  of  the  shell  and  the  propagation  medium  are 
used. 


4.1  Resultant  Pressure  Field  in  Front  of  Scatterer 


In  this  report,  the  resultant  pressure  magnitude  in  the  backscatter  direction  is 
computed.  This  pressure  field  (normalized  with  respect  to  the  inddent 
pressure  wave  amplitude)  may  be  writer  as 
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P  -  |P|/P^  -  |(P,/Po>*<P,-/Po)*(^/Po>l 

where  the  incident  and  elastic  scattered  pressures  dose  to  the  sphere  have  been 
determined  previously  to  be 

Pj/P<,  ■  (r,d)/PQ  -  expiiTrcosd) 


P.-/Po  "  ■  -  E  <2n>l)(t)"P„{cos^) 

n«o 


/n<^) 


^((T) 


fc„(Tr> 


(4^) 


P’^/P 


P^{r,t)/P„ 


^  (2n^l)(j)"P„(cosd)fc„(T»>) 


(4.4) 


where  depend  on  v,  ^  and 

t  ■  fcfl,  r  »  r/a,  v  =  h/a,  Q  =  t(c/Cp),  3^  =  Xn  =  n(n+l)  ^'*•5) 


> 


h„(x) 

/(pc)  —  z„ 

hUx) 


(  \ 


-«<P,Cp) 


¥ 


Tft^-(W3^)(v”^x^ 


(4.6) 


4.2  Resonance  Scattering  Theory 


The  resonance  scattering  theory  is  applied  in  order  to  analyse  the  backscattering 
cross-section  of  a  spherical  shell.  The  backscattering  cross-section  is  defined  as 
the  ratio  of  the  scattered  power  (referred  to  a  distance  of  1  metre),  to  the 
intensity  incident  on  a  unit  volume  [8],  and  nnay  be  written  as 


a 


■  4jt  I 


P.^(l,n) 


(4.7) 


Upon  normalizing  this  cross  section  with  respect  to  value  Jia^,  the  form  function 
f(n)  may  be  defined,  i.e. 


o 


2P„(a,n)  , 

I _ !! _ r 

Po 


\fM\^ 


where 


(4.8) 
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and  is  termed  the  n***  "partial  wave".  The  form  function  may  be  broken  up 
into  two  components;  the  backgrounds  and  the  resonances,  i.e. 


|/(JI)|  -  !/,»(«)  ♦/rU)  I  -  I  E  /ns-  ^  E  /n, 


(4.9) 


where  and  f„  are  the  modal  components  and 


2p,«(a,w)  .  -  (2n+l){-l)"j'(T) 

14,(7:)  I  -  iZ!!^ - |-|(2i)E. 


t)i„(t) 


(4.10) 


!/,(«) 


2P^(a,it) 


I  -  I  -(27pc)  5^ 


(2w*l)(-l)" 
n“o  (z„-.Z„)x"h'(x)^ 


(4.11) 


The  resonance  scattering  theory  is  useful  in  the  treatment  of  penetrable 
scatterers  in  that  it  establishes  connections  between  the  background  and 
resonance  features  observable  in  the  backscattering  cross-section  of  any  target. 


4,3  Target  Strength 


The  target  strength  in  the  far-held  may  be  defined  as 


TS  *  20  log|Q 


(4.12) 


J9-180’ 


20  log,o  l(r/P„)  1  Pjl.(r,6)  +  Pf  (r,0)  | 


where  P*'^  P*',  combine  to  give  the  scattered  pressure  in  the  far-field  and 

m 

(r/P^)p/,(r,d)  -  (fa/T)exp(iT7’)  ^  (2»t'*’I)P„(cosiJ)(/n(T)/ftn(T)] 

n"0 

(4.13) 


(r/P^)Pj(r,0)  -  -(i«/t)(pc)exp(iTf) 


(2n  +  l)P„(cos«) 


n-o  (z„+Z„)l(t))i^(t))2 


(4.14) 


4.4  Reflection  Factor 


Using  the  expressions  in  the  previous  section,  the  ratio  of  scattered  to  incident 
pressure  nnay  be  written  as 


T 


(4.15) 


The  1  /r  term  characterizes  the  spreading  of  the  scattered  wave,  whilst  the 
exponential  term  takes  into  account  the  phases  of  the  incident  and  scattered 
waves.  The  reflection  factor  can  be  made  independent  of  these  tenns  by 
defining  it  as 


to  yield 


“  (2b  ♦!)?  (cosi)) 

R  -  I  (i/t)  £  -  " - 

"•o  fc^(t) 


pc 


(t)2).'(t)(z„*2„) 


(4.17) 


The  reflection  factor  will  be  computed  as  a  function  of  angle  around  the  body, 
for  various  ka  values. 


4.5  Intensity  Field  Around  Scatterer 

Using  the  solution  to  the  resultant  pressure  held,  the  construction  of  two 
different  types  of  time-averaged  intensity  plots  are  possible,  including: 

(a)  Plots  of  the  intensity  magnitude  around  the  body  (both  "near-held"  and 
"far-held") 

(b)  Plots  of  the  intensity  vector  held  around  the  body  ("near-held"  only) 

For  both  type  (a)  and  (b)  plots,  the  intensities  for  the  scattered  and  resultant 
pressure  helds  are  computed.  The  radial  and  transverse  components  of  the 
intensity  vector  (averaged  over  one  p>eriod  of  the  wave)  are  dehned  as 

/,(r,«)  -  .1  Ar[P(r,d)H/(r,i9))  (4.18) 

^  PclP(r,d)iifl(r,6)l  (4.19) 
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where  Re  denotes  the  real  part,  •  the  complex  conjugate  and  u^,  the  velocity 
components  correspoitding  to  the  pressure  field  P  whether  it  lx;  a  scattered  or 
resultant  pressure.  The  magnitude  and  direction  of  the  intensity  vector  ate 
given  by  the  formulae 

^ +arctan(/^(r,d)//^(r,iJ)J  (4.20) 

t  -  arcunlz/il, 

where  X,  r  are  shown  in  Figure  1.  For  the  polar  intensity  plots  of  type  (a),  the 
intensity  value  is  first  divided  by  the  incident  intensity  =  P^^/^c.  The 
intensity  is  then  normalized  with  respect  to  the  maximum  value  that  it  obtains 
over  the  range  of  angles  from  0°  to  360“.  For  the  type  (b)  intensity  vector 
plots,  each  vector  has  first  been  divided  by  and  then  normalized  with  respect 
to  the  magnitude  of  the  largest  vector  in  the  plot.  The  square  root  of  the 
magnitude  is  then  calculated  and  multiplied  by  the  grid  spacing.  The  vectors 
are  plotted  over  the  cooidirute  range  -  2  S  (x/a)  S  2  and  0  s  (z/a)  i  2,  where  a 
is  the  mean  radius  of  the  spherical  scatterer. 

In  order  to  compute  the  tin\e-averaged  intensity,  the  velocity  field  associated 
with  the  pressure  field  must  be  found.  This  may  be  determined  from  Euler's 
equation,  vP  =  (i<op)u,  or  in  component  form 


1  -I  dP 
pcjlfciw 


1  -i  ap 
pc  kr  "So 


Using  the  "near-field”  and  "far-Seld"  pressure  expressions  given  in  sections  4.1 
and  4.2,  the  following  velocity  components  are  obtained. 


Velocity  due  to  incident  pressure  Pj 


(pc/P„)I«,],  -  ♦  cos6  exp(ifcrcos6) 
(pc/Pg)Iii^Jl  ■  -sini9exp(iiTC0sd) 


Velocity  due  tc  rigid  scattered  pressure 


Near-field: 


(pc/P„)lu,  ],.  -  £  an*V{if 


P„(cosi}) 


/n(^> 

- - hjxr) 
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*•  /  ^  {  T  ) 

(pc)/P^)(H^l,.- -(sine/t/')  (2n*l)(i)"’'P^(cosd)ii: _ h^{xf) 

h^(x} 

(425) 

Far-field: 


(pc/P^Xu/'],.  -  (i7tr)<'^^  5^  (2n*l)P„(cosO)  (426) 

"•o  ii^Ct) 


(pc/P^Xu^J,^  -  (-sind/(tf)^)e‘'^  Y  (2n *1  )P„^  (cosd)  .{T-JL-  (4-27) 


Velocity  due  to  the  radiated  pressure  Pj 
Near-field: 


(WP„Xu,^J 


E 


(2n  +  l)(i)"''P„(cosd)/i„^TO 
(2„*Z„X(T)A^t)P 


(4.28) 


(l/Po>l“el 


-(sintS/tr)  Y 

n«o 


(2n*l)(i)"*'P^(cose)fi„(tr) 

(2,-Z„)I(T)<(T))' 


(4.29) 


Far-field: 


(I/PoXm/J 
(l/P„)[u,^J  . 


-(-.7Tf)e"'  5: 

n»o 


(2n  ♦l)P„(cosd) 
(Z„»Zn)l(T)fi^(x)P 


(sind/(Tr)^)f‘^^  Y 
n»o 


(2n+l)Pn  (cosd) 
(z„.Z„)I(T)h'(t)]^ 


(430) 


(431) 


By  combining  equations  (2.1),  (2.17)  and  (3.16),  the  radial  and  ungential 
displacements  are  derived,  i.e. 
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(432) 


-  (2ii*l)(i)"P„(cosd) 

«»(«>  -  E  - - 


u(6) 


-  an(2n+l)(i)"{l-cos^«)’''^P'(cos«) 

- - 

(Zn 


(433) 


where 


,4 


o 


Po  -[p^(v^x^-l)*(l  ♦v)) 

(0(t)^'  l£i*-(l  ♦p^)(v*\,-l)I 


(434) 


Using  equations  (432)  and  (433),  the  deformed  shape  of  the  elastic  shell  may  be 
plotted. 


5.  Results  and  Discussion 


In  this  chapter,  plots  of  the  various  scattering  parameters  for  the  elastic 
spherical  shell  are  discussed  and  Gompared  with  those  of  the  rigid  sphere.  For 
some  plots,  only  four  specific  ka  values  will  be  considered,  namely  Jbi  =  1.14, 
1.44,  3.30  and  530.  These  values  correspond  to  the  n  =  2, 3  resonances,  as  well 
as  a  dip  and  a  maximum  in  foe  target  strength  sp>ectrum. 


5.1  Resultant  Pressure  Field 


Figures  2(a)  to  2(d)  show  the  normalized  resultant  pressure  as  a  function  of  r/a 
for  the  above  four  ka  values.  For  Figures  (a)  and  (b),  corresponding  to  the 
R  s  2,  3  resonances  of  the  shell,  radiation  scattering  far  outweighs  that  of  rigid 
body  scattering,  resulting  in  the  shell  pressure  being  many  times  that  of  its 
rigid  sphere  counterpart.  In  Figure  (c)  away  from  the  surface,  the  radiation 
scattering  aitd  rigid  scattering  interfere  destructively  thus  cancelling  each 
others'  effect,  whilst  for  Figure  (d)  the  shell  pressure  magnitude  is  comparable 
to  that  of  the  rigid  sphere  except  for  a  phase  difference  of  approximately  180°. 


Figure  2(a);  Resultant  Pressure  in  front  of  Elastic  Shell 


r/3 


Figure  2(b):  Resultant  Pressure  in  front  of  Elastic  Shell 


5.2  Resonant  Scattering  Theory 


As  mentioned  previously,  each  partial  wave  contained  in  the  form  function  of  a 
peitetrable  scatterer  can  be  tnoken  up  into  two  comporwnts,  a  smooth 
background,  and  a  set  of  modal  resonances.  Figures  3(a,b)  show  the  first  five 
modal  backgrounds  aiut  modal  resonances  for  the  thin  steel  spherical  shell. 

The  frequency  dependeiKe  of  the  form  function  contains  effects  which  arise 
from  four  distinct  types  of  contributions.  These  contributions  are  (1)  specular 
reflection,  (2)  creeping  waves,  (3)  symmetric  and  antisymmetric  Laird?  waves 
and  (4)  waves  bouiuting  within  the  inner  and  outer  boundaries  of  the  shell 
material. 

Specular  reflection,  which  is  effectively  the  rigid  body  reflection,  does  not 
produce  any  drastic  change  in  the  form  function  since  these  waves  do  not  travel 
along  any  part  of  the  surface.  Creeping  waves  are  diffracted  waves  thet 
originate  at  the  edge  of  the  shadow  boundary  and  reradiate  in  all  direc'.ons  as 
they  circumnavigate  the  body  in  the  fluid  medium.  Extrema  in  the  form 
function  due  to  waves  $catt«ed  within  the  elastic  boundaries,  can  only  occur 
for  high  frequencies. 


Background 

The  background  is  essentially  the  rigid  body  scattering  response,  i.e.  the 
response  of  the  scatterer  as  if  it  were  impenetrable,  and  is  due  to  creeping 
waves  travelling  around  the  shell.  The  background  is  the  contribution  due  to 
the  scatterers'  shape  only,  and  is  thus  independent  of  changes  in  shell  thickness 
and  material  composition. 


Modal  Resonances 

The  set  of  modal  resonances  is  tlw  radiation  scattering  response,  and  represents 
the  sound  radiation  from  the  shell  which  is  undergoing  forced  vibrations  due  to 
the  incident  wave  pressure.  This  energy  of  vibration  in  the  shell  goes  into  the 
creation  of  surface  (Lamb)  waves,  and  the  resonances  are  the  resultant  process 
of  phase  reinforcement  of  the  surface  waves  during  the  scattering  process  as 
th^  travel  around  the  shell. 

For  thin  shells  such  as  the  one  considered  in  this  study,  only  the  zeroth  order 
symmetric  Lamb  wave  is  dominant,  and  there  is  only  one  resonance  per  mode  n 
(see  Bg.  3(b)).  The  resulting  form  function  has  a  simple  structure  as  compared 
to  that  of  thick  shells,  where  many  Lamb  waves  interact  simultaneously, 
making  the  response  difficult  to  interpret.  Figure  4  shows  the  resultant 
backscattering  cross-section  for  the  steel  spherical  shell. 
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Figure  4;  Normalized  Backscattering  Cross-Section  of  a  Thin  Steel  Spherical  Shell. 


5.3  Target  Strength 

The  monostatic  tai]get  strength  spectre  of  the  thin  shell,  is  illustrated  in 
Figure  5,  where  the  dominant  resonant  peaks  are  labelled  by  their  individual 
modal  numbers.  The  resonant  peaks  of  this  shell  are  high,  with  the  n  &  2  peak 
being  almost  20  dB  above  the  ’rigid  body"  datum  plot.  The  oscillations  in  the 
T^d  sphere  curve  are  the  residt  of  interference  between  the  creeping  waves 
and  the  specular  reflectioa  There  are  no  resonances  for  either  ^e  n  =  0  or  the 
n  »  1  modes.  The  n  «  0  mode  which  is  termed  the  "breathing  mode” 
corresponds  to  a  uniform  erqsansion  and  contraction  of  the  spherical  shell, 
which  does  not  occur  for  pim  wave  ensonification.  The  n~  I  resonance 
occurs  at  zero  foequency  and  corresponds  to  an  urtdeformed  shell  moving  back 
and  forth  in  the  fluid  medium. 
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Figurt  5:  Monostatk  Target  Strength  of  an  Elastic  Shell. 


5,4  Reflection  Factor 

Figures  6(a}  to  6(d)  show  the  reflection  factor  for  the  four  ka  values  of  interest 
In  Figure  6(a),  the  scattered  field  is  synunetric  about  the  plane  containing  the 
shadow  boundary.  In  Figures  6(c)  and  6(d)  the  overall  appearance  of  the 
scattered  field  is  not  too  Cerent  from  that  for  a  rigid  sphm,  where  there  is  a 
strong  beam  of  interfering  radiation  behind  the  sphere.  There  are  however 
variations  in  the  magnitude  of  the  backscattered  echo  which  do  not  occur  with 
the  rigid  sphere  at  higher  frequencies. 
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Figure  6(tt):  Reflection  Factor  for  an  Elastic  Shell. 
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Figure  6<c):  Reflection  Factor  for  an  Elastic  Shell. 
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Figure  €<d):  Reflection  Factor  for  an  Elastic  Shell. 
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5.5  Intensity 


Figures  7  and  9  display  the  scattered  and  resultant  intensity  magiutudes  in  the 
near  held  as  a  function  of  airgle  around  the  scatterer.  The  scattered  response 
patterns  in  the  far-held  are  stown  in  Figure  8. 

The  scattered  and  resultant  inteisity  vector  fields  around  the  shell  are  shown 
in  Figures  10  and  11.  An  examination  of  the  resultant  fields  in  Figures  ll(a,b) 
clearly  shows  the  large  distortion  of  the  incident  sound  field  at  the  resoiuint 
frequencies.  At  ka  =  1.14  (n  s  2)  there  is  a  maximum  of  power  flow  into  the 
shell  at  the  poles,  and  a  iitaximum  of  power  flow  out  of  the  shell  at  the  equator. 
At  ka=  1.44  (n  =  3)  however,  there  is  a  mirumum  of  power  flow  into  the  shell 
at  the  poles,  but  still  a  maximum  of  power  flow  out  of  the  shell  at  the  equator. 

The  resultant  inteitsity  vector  fields  which  correspond  to  the  non-resonant  ka 
values  3.33  and  5.20  are  shown  in  Figures  ll(c,d).  Unlike  the  previous 
resonant  intensity  fields  which  display  complex  variations  in  pattern,  these 
fields  are  very  similar  to  those  generated  by  an  impenetrable  sphere.  The 
deformed  shapes  of  the  shell  at  the  four  ka  values  is  illustrated  in  Figures  12(a) 
to  12(d).  These  results  apply  to  deformations  which  occur  at  various  times 
during  the  excitation  period  and  demonstrate  the  types  of  results  which  can  be 
produced.  It  should  be  noted  that  the  deformations  as  shown  have  been  scaled 
by  the  factor  0.1  fl/{aKfi)lniax' 
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Figure  7(a);  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  7(b);  Scattered  Intensity  round  an  Elastic  Shell. 


Figure  7(c):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  7(d):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  8(a):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  8(b):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  8(c):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  8W):  Scattered  Intensity  around  an  Elastic  Shell. 
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Figure  9(a);  Resultant  Intensity  around  an  Elastic  Shell. 
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Figure  9(d):  Resultant  Intensity  around  an  Elastic  Shell- 
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j  Figure  10(a):  Scattered  Intensity  ll/Imax}^^  around  an  Elastic  Shell. 
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Figure  10(b):  Scattered  Intensity  ll/Imax]^^^  around  an  Elastic  Shell 
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Figure  10(e):  Scattered  Intensity  [I/Imax]^^  around  an  Elastic  Shell 
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Figure  10(d):  Scattered  Intensity  ll/Imax]^^  around  an  Elastic  Shell. 
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Figure  11(a):  Resultant  Intensity  Ullmaxl^^^  around  an  Elastic  Shell. 
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Figure  11(b);  Resultant  Intensity  around  an  Elastic  Shell. 
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Figure  life);  Resultant  Intensity  ll/Imaxp^  around  an  Elastic  Shell. 
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Figure  11(d):  Resultant  Intensify  lljlmoxfi^  around  an  Elastic  Shell. 
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Figure  12(a);  Deformed  shape  of  elastic  shell.  (Incident  pressure  Pg  =  10  Pa> 
(Scaling  factor  =  22  £406). 
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Figure  12(b):  Deformed  shape  dastk  shdl.  (Incident  pressure  s  10  Pa) 
(Scaling  factor  =  2.1  E+06) 


Figure  12(c):  Deformed  shape  of  dastic  shell.  (Inddent  pressure  Pg  =  10  Pa) 
(Scaling  factor  =  65  E+07). 
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Figure  12(d):  Deformed  shape  of  dastic  shell.  (Incident  pressure  Pg  =  10  Pa) 
(Scaling  factor  =  5.9  E+07). 


6.  Conclusion 


An  analysis  of  the  plane  wave  scattering  from  a  thin  steel  spherical  shell 
submerged  in  water  has  been  carried  out  and  results  are  presented  for  ha  values 
up  to  8.  Various  scattering  parameters  such  as  target  strength,  back-scattering 
ctoss-section,  reflection  factor  and  the  intensity  field  around  the  body  have  been 
computed  to  illustrate  the  solutions.  Results  from  a  rigid  sphere  of  comparable 
size  have  also  been  included  for  comparison  purposes.  At  values  of  ka 
corresponding  to  foe  shell  resonances,  intensity  vector  plots  show  large 
variations  in  foe  sound  field  around  foe  shell  and  target  strength  results 
indicate  echoes  of  up  to  20  dB  greater  than  that  obtained  from  a  similarly  sized 
rigid  sphere. 
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Resultant  pressure 
Incident  pressure  amplitude 
Incident  pressure 

Scattered  pressure  from  an  elastic  body 
Scattered  pressure  from  a  rigid  body  (infrnite  impedarKe) 
Radiated  pressure  from  an  elastic  b^y 
Pressure  at  the  shell  surface  (r  =  a) 

Backscattering  cross-section  of  target 

Form  function  in  backscatter  direction 

Modal  backgrounds  and  resonances  of  form  frmction 

Average  intensity  rtragrutude  of  the  incident  pressure 

Radial  and  tangential  components  of  the  intensity  vector 

Green's  function 

Angular  frequetrcy 

Normalized  frequency  of  irrcident  pressure  wave 

Natural  frequeiunes  of  vibration  for  a  spherical  shell  (f  =  I,  2) 

Acoustic  wavenumber 

Mean  radius  of  spherical  shell 

Speed  of  sound  in  water 

Speed  of  wave  in  a  shell 

Density  of  water 

Density  of  shell  material 

Young's  modulus  of  the  shell 

Poisson's  ratio  of  the  shell 

Shell  thickness 

Shell  thickness  factor 

Spherical  strains 

Spherical  stresses 

Spherical  coordinates 

Kinetic  energy 

Strain  energy 

Boundary  of  sphere  (r  =  a) 

Normal  velocity  component  of  resultant  pressure  at 
Normal  velocity  component  of  incident  pressure  at 
Normal  velocity  component  of  rigid  scattered  pressure  at  Sj, 
Normal  radiated  velocity  component  at  S„ 

Unit  outward  rtormal 
Modal  acoustic  impedance 
Modal  ntechanical  impedance 


52 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  UNCLASSIFIED 


CXXniMENT  CONTROL  DATA  SHEET 

REPORT  NO. 

MRL-TO-9M0 

ARNO. 

AR-006-354y 

REPORT  SECURiry  CLASSIHCATION 
Unclassified 

TITLE 

Sound  wave  scattering  from  an  elastic  spherical  shell 

AUTHOR® 

C.J.  Partridge 

CORPORATE  AUTHOR 

Materials  Research  Laboratory 

PO  Box  50 

Ascot  Vale  Victoria  3032 

REPORT  DATE 

January,  1993 

TASK  NO. 

NAV  89/020 

SPONSOR 

RAN 

FILE  NO. 

G6/4/8-4036 

REFERENCES 

8 

PAGES 

52 

CLASSIFICATION/UMITATION  REVIEW  DATE 

CLASSIFICATION/RELEASE  AUTHORITY 
Chief,  Underwater  Systems  Division 

SECONDARY  DISTRIBUTION 

Approved  for  public  release 

ANNOUNCEMENT 

Announcement  of  this  report  is  unlinuted 

KEYWORDS 

Wave  Scattering 
Acoustics 

Sound  Waves 

Spherical  Shell 

Pressure  Fields 

ABSTRACT 


The  Helmholtz  integral  equation  may  be  formulated  by  combining  the  scalar  wave  equation  and  Euler's 
equation  for  motion  within  a  fluid.  The  solution  to  this  integral  equation  yields  the  radiated  pressure 
from  a  submerged  vibrating  body.  The  solution  obtained  from  applying  Hamilton's  variational  principle 
governs  the  response  to  a  given  surface  pressure  excitation.  These  two  solutions  may  be  used  in 
conjunction  to  investigate  the  scatterii^  of  incident  sound  waves  from  shell-like  bodies.  In  this  report, 
the  scattering  from  a  submerged  steel  spherical  shell  having  a  thickness  to  radius  ratio  of  0.01  is 
investigated  for  ka  values  up  to  8. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
UNCLASSmED 


! 


Sound  Wave  Scattering  from  an  Elastic  Spherical  Shell 
C.).  Partridge 
(MRL-TK-91-10) 


DISTRIBUTION  LIST 


Director,  MRL 

Chief,  Maritime  Operations  Division 

Dr  Cl.  Sach 

C.J.  Partridge 

MRL  Information  Service 

Chief  Defence  Scientist  (for  CDS,  FASSP.  ASSCM)  (1  copv  onlv) 

Director,  Surveillance  Research  Laboratorv 

Director  (for  Library),  Aeronautical  Research  Laboratorv 

Director,  Electronics  Research  Laboratorv 

Head,  Information  Centre,  Defence  Intelligence  Organisation 

OIC  Technical  Reports  Centre,  Defence  Central  Library 

Officer  in  Charge,  Document  Exchange  Centre  (8  copies) 

Army  Scientific  Adviser,  Russell  Offices 

Air  Force  Scientific  Adviser,  Russell  Offices 

Navy  Scientific  Adviser,  Russell  Offices  -  data  sheet  only 

Scientific  Adviser,  Defence  Central 

Director-General  Force  Development  (Land) 

Senior  Librarian,  Main  Librarv  DST05 
Librarian,  MRL  Sydney 
Librarian,  H  Block 

UK/USA/CAN  ABCA  Armies  Standardisation  Rep.  c/-  DGAT  (8  copies) 
Librarian,  Australian  Defence  Force  Academy 
Counsellor,  Defence  Science,  Embassy  of  Australia  -  data  sheet  only 
Counsellor,  Defence  Science,  Australian  High  Commission  -  data  sheet  only 
Scientific  Adviser  to  DSTC,  C/-  Defence  Adviser  -  data  sheet  onlv 
Scientific  Adviser  to  MRDC,  C/-  Defence  Attache  -  data  sheet  onlv 
Head  of  Staff,  British  Defence  Research  and  Supply  Staff  (Australia) 

NASA  Senior  Scientific  Representative  in  Australia 

INSPEC:  Acquisitions  Section  Institution  of  Electrical  Engineers 

Head  Librarian,  Australian  .Nuclear  Science  and  Technology  Organi.sation 

Senior  Librarian,  Hargrave  Librarv,  Monash  University 

Library  -  Exchange  Desk,  National  Institute  of  Standards  and  Technology,  US 

Exchange  Section,  British  Library  Document  Supply  Centre 

Periodicals  Recording  Section,  Science  Reference  and  Information  Service,  UK 

Library,  Chemical  Abstracts  Reference  Service 

Engineering  Societies  Library,  US 

Documents  Librarian,  The  Center  for  Research  Libraries,  US 
CMDR  ).M.  Taubman,  Deputy  Director  Submarine  Warfare,  Canberra 


